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Quantum computers promise to solve certain 
problems exponentially faster than possible clas- 
sically but are challenging to build because of 
their increased susceptibility to errors. Remark- 
ably, however, it is possible to detect and cor- 
rect errors without destroying coherence by us- 
ing quantum error correcting codes 1 . The sim- 
plest of these are the three-qubit codes, which 
map a one-qubit state to an entangled three-qubit 
state and can correct any single phase-flip or bit- 
flip error of one of the three qubits, depending 
on the code used 2 . Here we demonstrate both 
codes in a superconducting circuit by encoding 
a quantum state as previously shown 3 ' 4 , induc- 
ing errors on all three qubits with some probabil- 
ity, and decoding the error syndrome by revers- 
ing the encoding process. This syndrome is then 
used as the input to a three-qubit gate which cor- 
rects the primary qubit if it was flipped. As the 
code can recover from a single error on any qubit, 
the fidelity of this process should decrease only 
quadratically with error probability. We imple- 
ment the correcting three-qubit gate, known as a 
conditional-conditional NOT (CCNot) or Toffoli 
gate, using an interaction with the third excited 
state of a single qubit, in 63 ns. We find 85±1% fi- 
delity to the expected classical action of this gate 
and 78 zb 1% fidelity to the ideal quantum pro- 
cess matrix. Using it, we perform a single pass of 
both quantum bit- and phase-flip error correction 
with 76 zb 0.5% process fidelity and demonstrate 
the predicted first-order insensitivity to errors. 
Concatenating these two codes and performing 
them on a nine-qubit device would correct arbi- 
trary single-qubit errors. When combined with 
recent advances in superconducting qubit coher- 
ence times 5 ' 6 , this may lead to scalable quantum 
technology. 

Quantum error correction relies on detecting the pres- 
ence of errors without gaining knowledge of the encoded 
quantum state. In the three-qubit code, the subspace 
of the two additional "ancilla" qubits uniquely encodes 
which of the the four possible single-qubit errors has oc- 
curred, including the possibility of no flip. Critically, er- 
rors consisting of finite rotations can also be corrected by 
projecting this syndrome, essentially forcing the system 
to "decide" if a full phase- or bit-flip occurred 2 . Previ- 
ous works implementing error correcting codes in liquid- 



7 9 and solid-state 10 NMR and with trapped ions 11,12 
have demonstrated two possible strategies for using the 
error syndromes. The first is to measure the ancillas 
and use a classical logic operation to correct the de- 
tected error. This "feed-forward" capability is challeng- 
ing in superconducting circuits as it requires a fast and 
high-fidelity quantum non-demolition measurement, but 
is likely a necessary component to achieve scalable fault- 
tolerance 2 ' 13 . The second strategy, as recently demon- 
strated with trapped ions 12 and used here, is to replace 
the classical logic with a quantum CCNot gate which 
performs the correction coherently, leaving the entropy 
associated with the error in the ancilla qubits. The CC- 
Not performs exactly the action that would follow the 
measurement in the first scheme: flipping the primary 
qubit if and only if the ancillas encode the associated 
error syndrome. 

The CCNot gate is also vital for a wide variety of ap- 
plications such as Shor's factoring algorithm 14 and has 
attracted significant experimental interest with recent 
implementations in linear optics 15 , trapped ions 16 , and 
superconducting circuits 17 ' 18 . Here we use the circuit 
quantum electrodynamics architecture 19 to couple four 
transmon qubits 20 to a single microwave cavity bus 21 , 
where each qubit transition frequency can be controlled 
on nanosecond timescales with individual flux bias lines 22 
and collectively measured by interrogating transmission 
through the cavity 23 . (The details of the device can be 
found in the Methods Summary and in Ref. 3.) Qubits 
are tuned to 6, 7, and 7.85 GHz, with the fourth at 
~ 13 GHz and unused (hereafter referred to as Q1-Q4). 
In this paper, we first demonstrate the three-qubit inter- 
action used in the gate, which is the logical extension of 
interactions used in previous two-qubit gates 3 ' 22 ' 24 , and 
demonstrate how this interaction can be used to create 
the desired CCNot. We then characterize its classical and 
quantum action and finally use the gate to demonstrate 
three-qubit error correction. 

Our three-qubit gate employs an interaction with the 
third excited state of one qubit. Specifically, it relies on 
the unique capability among computational states (<j z 
eigenstates) of | 111) (the notation \abc) refers to the ex- 
citation level of Q1-Q3, respectively) to interact with the 
non-computational state |003). As the direct interaction 
of these states is first-order prohibited, we first transfer 
the quantum amplitude of | 111) to the intermediate state 
1 102), which itself couples strongly to 1 003) . Calculated 
energy levels and time-domain data showing interaction 
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FIG. 1. Calculated energy spectra and time domain 
measurements of the interactions used in the three- 
qubit gate, (a) The energy spectrum of doubly excited states 
showing the avoided crossing between 1 01 1) and |002) (identi- 
cal to that between | 111) and 1 102) except for a 6 GHz offset) 
is shown with both (top) a numerical diagonalization of the 
system Hamiltonian and (bottom) a time-domain measure- 
ment as a function of the flux bias on Qi . (top) The frequen- 
cies for the involved eigenstates are blue and non-interacting 
eigenstates of similar energy are grey. The notation \abc) <g> \d) 
indicates the excitation level of each qubit and the cavity pho- 
ton number, respectively. When omitted, d = 0. (bottom) 
The state 1 01 1) is prepared and a square flux pulse of dura- 
tion t and amplitude V2 is applied. Coherent oscillations pro- 
duce a "chevron" pattern, with darker colors corresponding 
to population left in 1 002) . (b) The spectrum of triply excited 
states showing the avoided crossing between 1 102) and 1 003) 
as a function of the flux bias on Qi is characterized in the 
same way as above. 1 102) is prepared by first making | 111) 
and then performing the swap as described in Fig. 2. Many 
additional eigenstates are close in energy but are irrelevant 
because they do not interact with the populated states. A 
large avoided crossing between the relevant eigenstates that 
is used to produce an adiabatic three-qubit interaction hap- 
pens near 28 m$o- Extra lines near 31 m$o and 29 m$o are 
due to third-order interactions predicted by the Hamiltonian, 
as is the larger first-order interaction at 25 m$o, but their 
effect on the protocol in Fig. 2 is negligible. 



between 1 01 1) and 1 002) (which is identical to 1 1 1 1) and 
1 102) except for a 6 GHz offset) as a function of the flux 
bias on Q2 are shown in Fig. 1(a). Once the amplitude 
of I 111) is transferred to 1 102) with a sudden swap inter- 
action, three-qubit phase is acquired by moving Qi up in 
frequency adiabatically, near the avoided crossing with 
1 003). Figure 1(b) shows the avoided crossing between 
these states as a function of the flux bias on Q\. This 
crossing shifts the frequency of 1 102) relative to the sum 
of 1 100) and |002) to yield our three-qubit phase. The 
detailed procedure of the gate is shown in Fig. 2(a), tak- 
ing a total of 63 ns. Further details can be found in the 
Supplementary Information. 

We first demonstrate the gate by measuring its classical 
action. The controlled-controlled-phase (CCPhase) gate, 
which maps | 111) to — | 111), has no effect on pure com- 
putational states so we implement a CCNot gate by con- 
catenating pre- and post-rotations on Q2, as described 
in the Supplementary Information. Such a gate ideally 
swaps 1 101) and | 111) and does nothing to the remaining 
states. To verify this, we prepare the eight computational 
states, perform the gate, and measure its output with 
three-qubit state tomography 3 to generate the classical 
truth table. The intended state is reached with 85 ± 1% 
fidelity on average. This measurement is only sensitive 
to classical action, however, and a more thorough set of 
measurements is needed to fully characterize the gate. 

To complete our verification, we perform full quan- 
tum process tomography (QPT) on the CCPhase gate. 
In addition to detecting the action of the gate on quan- 
tum superpositions of computational states, QPT also 
detects non-unitary time evolution due to spurious cou- 
pling to the environment. It is done by preparing 64 
input states which span the computational Hilbert space 
and performing state tomography on the result of the 
gate's action on each state. As shown in Fig. 3, the fi- 
delity is found to be 78 ± 1% to a process in which the 
spurious two-qubit phase between Qi and Q3 is set to the 
measured value of 57 degrees (see the Supplementary In- 
formation for details on this phase and an explanation of 
why it is irrelevant here). Due to this extraneous phase, 
the phase gate is most accurately described as a CC-e^Z 
gate (Z is a Pauli operator 2 ). The infidelity is consistent 
with the expected energy relaxation of the three qubits 
during the 85 ns measurement, with some remaining er- 
ror owing to qubit transition frequency drift during the 
90 minutes it takes to collect the full dataset. 

With our CCPhase gate in hand, we now demonstrate 
three-qubit error correction. Both the phase- and bit-flip 
codes begin by encoding the quantum state to be pro- 
tected in a three-qubit entangled state 2 by using condi- 
tional phase (CPhase) gates, as shown in Fig. 4(a). The 
two codes differ only by single-qubit gates applied after 
entanglement in the encoding step. For quantum states 
on the equator of the Bloch sphere, the resulting encod- 
ing is a maximally entangled three-qubit GHZ state 3 ' 4,26 
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FIG. 2. Three-qubit gate pulse sequence and classi- 
cal action, (a) The frequency of the three qubits during the 
gate as a function of time. First, Qi is moved suddenly into 
resonance with the avoided crossing shown in Fig. 1(a) to 
coherently transfer the population of 1 1 1 1 ) to 1 102) (and also 
1 011) to 1 002)) in 7 ns. Fine adjustments in the first point of 
the pulse compensates for finite pulse rise time and temporal 
precision. Qi is then moved suddenly further up in frequency, 
to where its two-qubit phase with Q3 is cancelled during the 
gate by accumulating a multiple of 2tt. Qi is then moved 
up adiabatically to initiate the interaction between 1 102) and 
1 003). The duration and amplitude of this pulse is tuned to 
acquire a three-qubit phase of exactly tt. The population in 
1 102) is then transferred back to 1 1 1 1) by reversing the swap 
procedure. Finally, the two-qubit phase between Q\ and Qi 
is cancelled with an additional adiabatic interaction, which is 
sped up with a tt pulse on Qi at 37 ns. Here, this tt pulse 
is explicitly undone after the gate, but when it is used for 
error correction the following pulse is compiled together with 
other post-rotations. The two-qubit phase between Qi and 
Q3 is uncontrolled, making this a CC-e^Z gate, (b) A CC- 
Not gate is made by appending to the phase gate pre- and 
post-rotations on Qi as described in the Supplementary In- 
formation. Its classical action is measured by preparing the 
eight computational basis states and performing state tomog- 
raphy on the result of applying the gate to them. The projec- 
tion of these measurements with the computational basis is 
taken to generate the truth table and is plotted. The fidelity 
to the expected action, where only the states 1 101) and 1 1 1 1) 
are swapped, is 85 =b 1%. 



which we independently measure to have a state fidelity 
of 89 ± 1%. Once the state is encoded, a single error 
of the chosen type on any of the qubits can be detected 
and corrected. The error syndrome is decoded by revers- 
ing the encoding sequence, leaving the ancilla qubits (Qi 
and Qs) in a state indicating which error occurred. For 
a full flip, they will be in a computational state. In par- 
ticular, both ancillas will be excited if the primary qubit 
(Q2) was flipped, and so the application of the CCNot 
gate will correct it. As detailed in the Supplementary 
Information for the case of bit-flip errors, an arbitrary 
rotation on any single qubit about the protected axis can 
be also encoded, detected, and reversed. 

In real physical systems, errors will occur at approx- 
imately the same rate on all constituent qubits. The 
correction scheme will succeed, therefore, when the sys- 
tem projects to zero or one errors. The probability of 
more than one error occurring is 3p 2 — 2p 3 , where p is 
the single-qubit error rate 2 , and so the fidelity of error 
correction should be 1 — 3p 2 +2p 3 . For a scheme with gate 
fidelity limited by decoherence, these coefficients will be 
smaller, but crucially, any linear dependence on p will 
be strongly suppressed. As shown in Fig. 4(b), we mea- 
sure the process fidelity of the phase-flip error correction 
scheme as a function of p by encoding four states which 
span the single-qubit Hilbert space and performing state 
tomography on the procedure's output. We compare this 
to the case of no error correction in which identical single- 
qubit rotations are applied to Q2 but the ancillas are not 
involved (and with appropriate delays to have the same 
total procedure duration). Whereas without error cor- 
rection we find a purely linear dependence on p, with 
the correction applied the data is extremely well mod- 
eled by only quadratic and cubic terms, demonstrating 
the desired first-order insensitivity to errors. 

We have realized both bit- and phase-flip error cor- 
rection in a superconducting circuit. In doing so, we 
have tested both major conceptual components of the 
nine-qubit Shor code 1 , which can protect from arbitrary 
single-qubit errors by concatenating the bit- and phase- 
flip codes. The implementation relies on our efficient 
three-qubit gate which employs non-computational states 
in the third excitation manifold of our system, demon- 
strating that the simple Hamiltonian of the system accu- 
rately predicts the dynamics even at these high excita- 
tion levels. The gate takes approximately half the time 
of an equivalent construction with one- and two-qubit 
gates. We expect it to work between any three nearest- 
neighbor qubits in frequency regardless of the number 
of qubits sharing the bus, as interactions involving other 
qubits will be first-order prohibited. 

We thank G. Kirchmair, M. Mirrahimi, I. Chuang, 
and M. Devoret for helpful discussions. We acknowl- 
edge support from LPS/NSA under ARO Contract No. 
W911NF-09-1-0514 and from the NSF under Grants No. 
DMR-0653377 and No. DMR-1004406. Additional sup- 
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FIG. 3. Quantum process tomography of the three-qubit phase gate. Absolute values of the elements of the (a) ideal 
and (b) measured process matrices. Data is collected by preparing 64 input states which span the three-qubit Hilbert space, 
applying the phase gate to them, and measuring the resulting density matrix with state tomography. The process matrix x 

of the operator O is related to these data by p ou t — O(p-m) = J2m n =i XmnAmpinA^, where pi n is the density matrix of the 
input state, p ou t is the measured output, Ai is an operator basis spanning the three-qubit operator space, here chosen to be 
the tensor products of three Pauli matrices, and N = 3 qubits 2 . The operator basis is ordered as in Ref. 3 and is explicitly 
written in the Supplementary Information. The ideal nonzero bars along the left edge are III, IIZ, IZI, ZII, IZZ, ZIZ, ZZI, and 
ZZZ. The fidelity of the operation / = Tr[x e x P tXthry] = 78 =b 1%. The ideal process matrix is calculated with the uncorrected 
phase between Qi and Q3 set to its measured value of 57 degrees, which is irrelevant for our implementation of quantum error 
correction because the ancilla qubits would be reset to their ground state 25 for a repeated cycle of correction. The fidelity to 
the "true" CCPhase gate, where the Q1-Q3 phase is set to 0, is 69 =b 1%. 



port provided by CNR-Istituto di Cibernetica, Pozzuoli, 
Italy (LF) and the Swiss NSF (SEN). 

METHODS 
Hamiltonian parameters 

The Tavis-Cummings Hamiltonian describing our sys- 
tem with four transmon qubits is 

H — huj c a^a + 

q=l j=0 j,k=0 

Here, H is Planck's reduced constant, u c is the bare cav- 
ity frequency, is the transition frequency for trans- 
mon q from ground to excited state j, and = g q rijk, 
with g q a bare qubit-cavity coupling and njk a cou- 
pling matrix element, uff and rijk depend on trans- 
mon charging (Ec q ) and Josephson (Ej q ) energies 27 . 
Flux dependence comes from E Jq = i^ ax |cos(7r<I>g/<l>o)|, 
with $> q the flux through the transmon SQUID loop 
and $0 is the flux quantum. A linear flux-voltage re- 
lation <fr q — Y^t=i a qi^ , i + ^g,o describes crosstalk and 
offsets. Spectroscopy and transmission data as a func- 
tion of flux bias gives uj c /2tt = 9.070 GHz, Ef™/h = 
{33, 35, 26, 57} GHz (from Q 1 to Q 4 ), g q /2n « 220 MHz, 



and Ec q /h ~ 330 MHz. The measured qubit lifetimes 
for Q1-Q3 are T± = (1.3,0.9,0.7) /is and coherence times 
T 2 * = (0.5,0.6, 1.3) /is respectively. 



Qubit rotations and gate calibration 

Arbitrary qubit rotations around the x- and ?/-axis of 
the Bloch sphere are performed with pulse-shaped reso- 
nant microwave tones. Rotations around the z-axis are 
done by rotating the reference phase of subsequent x and 
y pulses. One-qubit dynamical phases resulting from flux 
excursions are measured with modified Ramsey exper- 
iments comparing the phase difference between an un- 
modified prepared state and that same state after a flux 
pulse and are cancelled with z rotations. Two- and three- 
qubit phases are measured with a similar Ramsey exper- 
iment comparing the phase difference acquired when a 
control qubit is in its ground and excited state. For ex- 
ample, the two-qubit phase between Q 2 and Q3 is mea- 
sured by preparing Q 3 along the 2/- axis and Q2 either in 
its ground or excited state and then performing the flux 
pulse in both cases. The single-qubit phase of Q3 is the 
same for both states, and so the two-qubit phase is di- 
rectly measurable as their phase difference. All phases 
are initially tuned to within one degree, limited by the 
resolution of control equipment and drifts of system pa- 
rameters such as the qubit transition frequencies. 
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FIG. 4. Three-qubit phase-flip error correction 
scheme and demonstration of first-order insensitivity 
to errors, (a) The error correction protocol starts by entan- 
gling the two ancilla qubits with the primary qubit through 
the use of two CPhase gates (vertical lines terminating in solid 
circles). The number adjacent to each indicates which state 
receives a phase shift. A tt/2 rotation on the primary qubit 
is then performed, making this a phase- flip error correction 
code. If we wished to protect from bit flips, the two ancilla 
qubits would instead be rotated 2 . We perform errors on all 
three qubits simultaneously with z-gates of known rotation 
angle, which is equivalent to phase- flip errors with probability 
p = sin 2 (0/2). The encoding is then reversed, leaving the an- 
cillas in a state indicating which single-qubit error occurred. If 
an error has occurred on the primary qubit, the CCNot gate 
implemented with our CCPhase gate (represented by three 
solid circles linked by a vertical line) at the end of the code 
will reverse it. We then perform three-qubit state tomogra- 
phy to measure the result, (b) The fidelity of the protected 
qubit process matrix to the identity operation is plotted as 
a function of p. As the code corrects only single-qubit er- 
rors, it will fail if more than one error occurs, which happens 
with probability 3p 2 — 2p 3 . These coefficients are reduced for 
processes with finite fidelity The process fidelity is fit with 
/ = (0.76 ± 0.005) - (1.46 ± 0.03)p 2 + (0.72 ± 0.03)p 3 . If a 
linear term is allowed, its best-fit coefficient is 0.03=1=0.06. We 
compare this to the case of no error correction to simulate the 
improvement, (insets) The constituent state fidelities of the 
four basis states used to produce the process fidelity data for 
the case of (right) error correction and (left) no correction. 
The state \+Y) is immune to errors because its encoded state 
is an eigenvector of two-qubit phase flips. 
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DETAILS OF THREE-QUBIT PHASE GATE 

To understand the physical mechanism behind our 
three-qubit gate, it is useful to first review how two- 
qubit CPhase gates are commonly implemented in the 
cQED architecture. The avoided crossing between the 
first excited state of two transmon qubits (|11)) and the 
second excited state of one (|02)) can be employed both 
adiabatically 1 or suddenly 2 ' 3 to produce a CPhase gate. 
If the system is adiabatically tuned to the vicinity of 
the avoided crossing, the system will remain fully in the 
eigenstate that maps to the computational subspace away 
from the crossing. However, the interaction causes the 
energy level of the |11) state to differ from the sum of 
the 1 01 ) and 1 10) states, advancing the quantum phase 
of 1 11) relative to the quantum phases of 1 10) and |01), 
and entangling the qubits 1 . Alternatively, if the avoided 
crossing is small enough, the gate can be executed sud- 
denly. At the avoided crossing, the eigenbasis consists of 
the symmetric (|+)) and antisymmetric (|— )) superpo- 
sitions of 1 11) and |02). Starting from 1 1 1 ) and moving 
suddenly, the wavefunction will no longer be in an en- 
ergy eigenstate, but rather an equal superposition of |+) 
and |— ) whose relative phase will advance with the inter- 
action strength, oscillating between the undressed states 
1 11) and 1 02). Waiting one full period will return the state 
to the computational basis, but with a phase difference 
of 7r, at which point it can be moved suddenly away from 
the avoided crossing to return to the undressed eigenba- 
sis. 

Any two-qubit phase gate can be described in terms 
of one- and two-qubit phases, which is a helpful ab- 
straction to understand the different sources of phase 
delay. In this language, a two-qubit phase gate maps 
|00) |00), |10) -> e^ 10 |10), |01) -> e^ 01 |01), and 
|11) e^io+^oi+^iOin^ Here </> 10 and 0oi are single- 
qubit phases given by the time-integrated detuning of a 
qubit transition frequency from its nominal bias point, 
and 0n is a two-qubit phase which can be generated as 
described above. Note that the state 1 1 1 ) suffers phases 
from all of these sources, and so measuring the two-qubit 
interaction can only be done with separate measurements 
of 0oi an d 0io- Combining an interaction generating 
0ii = 7r with appropriate one-qubit phases can yield a 
CPhase gate conditioned on any target computational 
state 1 . 



The three-qubit phase of our C CPhase gate arises from 
an interaction between | 111) and 1 003) , in analogy to the 
two-qubit case. In the same way that the CPhase gate 
requires two excited qubits to access |02), the CCPhase 
gate requires three excited qubits to access |003). The 
direct interaction between these states is first-order pro- 
hibited because the states differ by a change of two exci- 
tations, so we achieve an effective interaction by using a 
state which couples strongly to both: 1 102) . As described 
in the main text, the gate is initiated with a full coherent 
transfer of the population of 1 1 1 1) to 1 102) (and 1 01 1) to 
1 002)) by suddenly approaching the crossing and oscillat- 
ing between the two dressed eigenstates for exactly half 
the splitting period before moving the qubit suddenly 
further up in frequency. Small errors due to imperfect 
timing and finite pulse rise time are corrected by shaping 
the qubit's trajectory to minimize residual |002) popula- 
tion after swapping into the state and back. Thanks to 
the strong coupling (67 MHz) of these states, the trans- 
fer takes only 7 ns. Once the amplitude of the target 
state is transferred to 1 102) , the three-qubit phase is ac- 
quired by moving Qi up in frequency, near the avoided 
crossing with |003). The interaction strength between 
these states is large (121 MHz) so we choose to acquire 
this three-qubit phase adiabatically. As in the case of 
two qubits described above, 1 102) experiences a frequency 
shift relative to its constituents ( 1 100) plus |002)) due to 
the avoided crossing, yielding the three-qubit phase. 

Again extending the two-qubit case, a three-qubit 
phase gate can be parametrized with seven unique 
phases. Three are one-qubit phases (0ooi, 0oio, and 
0ioo), three are two-qubit phases (0on, 0ioi? an d 0no)> 
and one is a three-qubit phase (0m). Ideally, our gate 
procedure would provide full control over them all, but 
we make a simplifying assumption based on the intended 
application of the gate. During error correction, errors 
are transferred from the protected qubit to the ancil- 
las, which are then reset either through measurement- 
conditioned pulses or by coupling them to a dissipative 
bath 5 . As their final state does not matter, we are free to 
choose any one two-qubit phase to remain uncorrected: 
here, the non- nearest neighbor interaction given by 0ioi, 
which is most challenging to control. This implies that 
Q2 will be the target of our error correction scheme, with 
Qi and Q3 acting as the ancillas. The remaining two- 
qubit phases 0on and 0no, however, must be set to zero. 
The former can be easily corrected via fine-tuning of the 
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Pauli operator Pauli operator 

FIG. SI. Reconstructed density matrices of the result of applying the CCNot gate to computational states. 

Each computational state is prepared and the CCNot gate described below is acted on it. Ideally, all states would have 
nothing happen to them except for 1 101) and 1 1 1 1) , which swap. The computational state prepared and what it should map 
to is indicated on each tomogram, which is visualized as the Pauli set as found in Ref. 2. The Pauli set of the ideal state is 
superposed (open bars) and is reached with fidelity 95.4, 95.2, 87.1, 85.2, 84.4, 82.2, 78.1, and 75.7% for each of the states, 
respectively. Note that the lack of spurious two- and three-qubit correlations indicates that there is no significant loss of 
population from the computational subspace. 



frequency of Q2 while in the \x02) state {x = 0, 1), where 
it acquires that phase very rapidly and so can be made to 
be an integer multiple of 2tt. The latter phase, however, 
must be explicitly corrected with an additional adiabatic 
phase gate. Because this angle is small, it is advanta- 
geous to 7r-pulse Q2 prior to the adiabatic interaction, 
reversing the direction of phase evolution and reducing 
the overall correction time. The action of the gate is to 
set 0ioo, 0oio, 0ooi, 0iio, and O n to zero, m to 7r, 
with 0ioi measured to be 57 degrees. 

A CCNot gate is constructed by appending a tt/2 and 
a — 7r/2 pulse on Q 2 before and after the phase gate de- 
scribed in Fig. 2(a) of the main text. For the two input 
states where both ancillas (Qi and Q3) are excited, there 
is an effective tt phase shift applied to the second pulse, 
and so the two pulses add together to a full rotation. 
Other states are not shifted so the pulses cancel. The 
full state tomograms of the output of CCNot with the 
eight computational states as input are shown in Fig. SI. 
These tomograms are used to derive the classical truth 
table shown in Fig. 2(b) of the main text. 

All relevant phases are controlled to one degree or bet- 
ter, with their accuracy set by the voltage resolution of 
our arbitrary waveform generator. This implies a "quan- 
tum phase fidelity" as defined in a recent work 6 in excess 
of 99% to the six relevant phases. This metric is not 



very sensitive to even phase errors, however, and so true 
quantum process fidelity should be used whenever pos- 
sible. For example, the phase fidelity of an ideal (e.g. 
decoherence-free) CCPhase gate with a maximal single 
qubit phase error of 0ioo = 7r to a CCPhase in which 
0ioo = is 85.7% despite the fact that the quantum pro- 
cess fidelity between those gates reveals them to be nearly 
orthogonal, with a fidelity of 6.3%. The real and imagi- 
nary parts of the measured process matrix are shown in 
Fig. S2. 

BIT-FLIP ERROR CORRECTION 

In order to illuminate the error projection process, we 
demonstrate bit-flip error correction with errors on only 
one qubit at a time. As in the phase- flip case described 
in the main text, bit-flip correction begins by encoding 
a quantum state in a three-qubit entangled state 7 using 
two sudden CPhase gates paired with appropriate single- 
qubit rotations 2 . The difference between the phase- and 
bit-flip codes lies only in the rotations performed after 
the entanglement. Instead of rotating the primary qubit 
as in phase-flip correction, the ancillas are tt/2 pulsed. 
For simplicity, here we only measure state tomography 
(rather than process tomography) of the state most sen- 
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FIG. S2. Real and imaginary parts of ideal and measured process tomography % matrix. These data were collected 
as described in Fig. 3 of the main text. There only the absolute value is shown, but here the full real and imaginary parts 
are reproduced. The order of operators here and in Fig 3. of the main text is as follows: III, IIX, IIY, IIZ, IXI, IYI, IZI, XII, 
YII, ZII, IXX, IYX, IZX, IXY, IYY, IZY, IXZ, IYZ, IZZ, XIX, YIX, ZIX, XIY, YIY, ZIY, XIZ, YIZ, ZIZ, XXI, YXI, ZXI, 
XYI, YYI, ZYI, XZI, YZI, ZZI, XXX, YXX, ZXX, XYX, YYX, ZYX, XZX, YZX, ZZX, XXY, YXY, ZXY, XYY, YYY, ZYY, 
XZY, YZY, ZZY, XXZ, YXZ, ZXZ, XYZ, YYZ, ZYZ, XZZ, YZZ, and ZZZ. We do not make use of the maximum-likelihood 
estimator commonly used to require the physicality of x matrix so that the reported elements of x an d the fidelity are linearly 
related to the raw measurements 4 . The uncertainty of the fidelities reported in the main text is given by the standard deviation 
of six repeated measurements of the full process matrix. 



sitive to bit-flip errors of the chosen type: the positive 
eigenstate of <j x . The state is now in a protected subspace 
which can recover from any single spurious y rotation on 
any of the three qubits. 

We perform intentional rotations on one of the qubits 
with a varying rotation angle 0. In the framework of 
error correction 7 , errors are decomposed as probabilis- 
tic full bit flips rather than continuous rotations, and 
so our partial rotations can be instead seen as varying 
the probability of a full bit flip. In the case of the an- 



cillas, this ^/-rotation would normally be placed between 
a positive and negative ir/2 x-rotation associated with 
turning a CPhase into a CNot, and so we compile these 
three single-qubit gates into one z-gate. After the error 
has occurred, we disentangle the three-qubit state, effec- 
tively encoding an error code in the ancillas. This code 
will leave both ancillas excited if and only if a bit flip has 
occurred on the primary qubit, which is then reversed by 
our CCNot gate. At this point, the entropy of the error 
is stored in the ancillas, which should be reset via cou- 
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FIG. S3. Bit-flip error correction, (a) Bit-flip error cor- 
rection gate sequence. This differs from Fig. 4(a) of the 
main text by only single-qubit rotations. Note also that for 
simplicity here we do only state tomography on the output 
of the process on the maximally affected state, so the state 
preparation is as indicated for Q2. (b) State fidelity to the 
created state after performing a single error on only one of 
the qubits, with and without error correction. Ideally, the 
curves would be flat lines at unit fidelity. Finite excited-state 
lifetimes cause oscillations and displacement down as the er- 
rors change the excitation level of the system, (c) Two-qubit 
density matrices of the ancillas after each of the four possible 
errors has occurred. The fidelity of each of these states to the 
ideal error syndromes are (81.3%, 69.7%, 73.1%, 61.2%). 



pling to a cold bath 5 if we were to loop the code. In Fig. 
S3(b), we plot the fidelity of the protected qubit to its 
prepared state after the error correction procedure has 
been applied to errors on all three qubits and also if no 
error correction is done. Ideally, these curves would be 
flat and with unit fidelity, but because of qubit decay and 
the varying excitation level of the qubits depending on 
the error performed, the curves show a small oscillation 
centered at 75.7% fidelity. 

We also measure the density matrices of the ancilla 
qubits after the four possible full bit-flip errors (no er- 
ror and bit flips on each of the three qubits). The 
code should ideally encode each of those four possible 
errors as one of the four computational states of the two- 
qubit ancilla subspace. As shown in Fig. S3(c), the 
measured density matrices of the ancillas do indeed en- 
code the error syndrome as expected, albeit with finite 
fidelity. The measured state fidelity to the ideal syn- 
drome is / sy ndrome = (81.3%, 69.7%, 73.1%, 61.2%) to the 
states ( 1 00) , |10), |01), |11)) encoding no error and errors 
on Q3, Qi, and Q2 respectively. In the case of a finite 
rotation, the ancillas will instead be in a superposition 
of error and no error, and so the gate will coherently 
correct the primary qubit, acting only on the subspace 
where the implicated error occurred. This action is free 
to be done by measurement and conditional feed-forward 
pulses (that is, by projecting the qubits with measure- 
ments), however, and indeed simplifies the requirements 
for fault-tolerance if it is done that way 7 . 
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